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In  the  smart  grid  environment,  the  flexible  and  diverse  distributed  generation  (DG)  and  microgrid  (MG) 
are  attracting  considerable  attention.  There  are  many  key  management  and  optimization  issues  involved 
in  smart  grid.  As  an  important  part  of  smart  grid  optimizations,  the  optimal  load  distribution  of  MG 
contributes  to  the  efficient  operation  of  MG  in  the  smart  grid  environment.  However,  traditional  optimal 
load  distribution  models  of  large-scale  power  generation  systems  are  not  fully  applicable  to  MG  for  their 
distinct  characteristics.  In  this  paper,  we  first  introduce  the  MG  in  smart  grid  and  analyze  its 
characteristics.  Then,  we  present  a  review  of  the  optimal  load  distribution  models  of  MG  in  the  smart 
grid  environment,  and  point  out  the  deficiencies  of  the  existing  models.  Finally,  a  comprehensive 
optimal  load  distribution  model  of  MG  both  in  objective  functions  and  constraints  is  established  and 
discussed. 
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1.  Introduction 

With  the  rapid  development  of  power  generation  technologies 
and  the  increasing  demand  for  electricity,  large-scale  centralized 
*  Corresponding  author  power  generation  systems  were  widely  implemented,  in  which 
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the  disadvantages  of  large-scale  centralized  thermal  power  gen¬ 
eration,  such  as  disadvantages  in  the  aspects  of  economic  effi¬ 
ciency,  environmental  protection,  security  and  reliability,  and 
control  agility,  were  gradually  emerging.  People  began  to  re¬ 
evaluate  the  large-scale  centralized  power  generation  methods, 
and  the  idea  of  small-scale  distributed  generation  was  presented 
[1],  After  the  late  1990s,  increasing  research  efforts  have  focused 
on  the  distributed  generation  (DG)  and  its  economic  benefits  [2- 
5],  DG,  also  known  as  dispersed  generation  or  distributed  energy 
supply,  is  a  new  energy  production  and  supply  model  [6-8].  In  DG, 
relatively  small  power  generation  units  are  arranged  near  the  user 
load  sites.  Therefore,  compared  with  the  large-scale  centralized 
power  generation,  the  main  advantages  of  DG  are  its  dispersion  in 
locations  and  flexibility  in  power  supply.  In  addition,  DG  can  also 
improve  the  reliability  of  power  supply  system  when  collaborating 
with  large-scale  power  generation  systems  [9],  However,  DG  also 
has  some  disadvantages.  The  main  power  grid  always  employs 
methods  like  restriction  and  isolation  to  reduce  the  impact  from 
DG.  With  the  development  of  smart  grid,  the  concept  of  microgrid 
(MG)  [10-12]  was  proposed  in  order  to  coordinate  the  contra¬ 
dictions  between  DG  and  the  main  power  grid.  MG  is  an  important 
power  generation  mode  in  smart  grid  [13-15], 

Currently,  MG  is  primarily  in  the  experimental  research  stage 
[16-18],  There  are  still  many  obstacles  in  terms  of  technology, 
control  and  management  that  needs  to  be  overcome  before  MG 
can  be  widely  applied  [19-21].  The  economic  optimization  of  MG 
is  one  of  the  key  management  issues  which  can  affect  its  operating 
efficiency  and  user  acceptance  [22-24],  Three  aspects  of  MG 
economics  were  summarized  in  [25].  Solar  power,  wind  power 
and  other  renewable  energy  power  generation  methods  are 
widely  adopted  in  MG,  which  makes  its  environmental  benefits 
apparent  [26].  Also,  power  generated  in  MG  can  be  supplied  in  a 
safer,  more  reliable,  and  more  personalized  way  by  means  of 
electricity  users  division  and  load  classification  [27,28].  Currently, 
there  have  been  some  research  efforts  in  the  economic  optimiza¬ 
tion  of  MG  [29-32], 

The  various  optimization  issues  are  the  key  parts  of  MG 
economics,  such  as  the  location  optimization  [33-35]  and  capacity 
optimization  [36-38]  of  distributed  generators  in  order  to  deter¬ 
mine  the  location  and  capacity  of  each  generator.  The  optimal  load 
distribution  is  an  important  optimization  problem  to  support  the 
efficient  operation  of  MG  [39],  The  optimal  load  distribution  of  MG 
is  to  achieve  multi-objectives,  including  low  operating  cost,  low 
emissions,  high  reliability,  high  power  quality,  and  low  line  loss, 
etc.,  while  meeting  various  system  operating  constraints.  In  this 
paper,  we  will  review  some  existing  optimal  load  distribution 
models  of  MG,  and  point  out  the  deficiencies  of  the  existing 
models.  Then  a  general  and  comprehensive  model  of  MG  optimal 
load  distribution  is  proposed  and  discussed. 

The  remainder  of  this  paper  is  structured  as  follows.  A  brief 
introduction  of  smart  grid  and  MG,  as  well  as  the  optimal  load 
distribution  of  MG  are  described  in  Section  2.  Then,  the  optimiza¬ 
tion  mathematical  models  of  MG  optimal  load  distribution,  both 
single-objective  and  multi-objective,  are  reviewed  in  Section  3.  In 
Section  4,  we  propose  a  comprehensive  model  of  MG  optimal  load 
distribution,  including  different  objective  functions  and  various 
constraints.  Finally,  conclusions  are  drawn  in  Section  5. 


2.  Optimal  load  distribution  of  MG  in  the  smart  grid 
environment 

2.J.  MG  in  the  smart  grid 

To  achieve  the  goal  that  the  operation  of  power  systems 
becomes  safer,  more  reliable,  more  environmentally  friendly,  more 


flexible,  more  controllable,  and  more  cost-effective,  the  United 
States  and  some  European  countries  proposed  the  concept  of 
“smart  grid”  [40-42],  Currently,  many  countries  have  put  smart 
grid  as  one  of  the  national  strategies,  and  considerable  research 
work  have  been  carried  out  [43-45],  Smart  grid  [46-48]  can  be 
generally  considered  as  an  intelligent  power  grid  system  which 
integrates  the  energy  flow  and  information  flow  by  advanced 
information  technology,  sensor  technology,  automatic  control 
technology  and  scientific  management  methods,  etc. 

MG  is  an  important  part  of  smart  grid  implementation  and 
application.  A  variety  of  distributed  generators  and  energy  storage 
devices  can  assess  to  the  main  power  grid  in  the  smart  grid 
environment.  They  can  form  a  “virtual  power  plant”  through  the 
interconnection  on  all  voltage  levels  by  electrical  communication 
systems.  Many  researchers  studied  the  roles  of  distributed  gen¬ 
erators  and  energy  storage  devices  in  smart  grid.  Zhang  [49] 
described  a  framework  for  the  operation  and  control  of  smart  grid 
with  DG  and  Flexible  AC  Transmissions  (FACTS),  and  proposed  a 
global  coordinated  strategy  for  voltage  control.  Mohd  et  al.  [50] 
pointed  out  that  distributed  energy  storage  systems  with 
advanced  power  electronics  can  play  an  important  role  in  power 
supply  systems  and  lead  to  many  financial  benefits.  They  also 
studied  the  topologies,  control  and  flexibility  of  the  main  power 
grid  with  energy  storage  systems.  Chowdhury  et  al.  [51]  studied 
the  operation  and  control  of  distributed  generators  in  power 
islands  by  simulation,  and  an  appropriate  bus  controller  was 
designed.  They  also  analyzed  the  flexibility,  safety  and  power 
quality  of  DG.  These  research  efforts  mainly  focused  on  the 
technical  and  control  perspective  of  DG  operation  in  the  smart 
grid  environment. 

2.2.  The  characteristics  of  MG  optimal  load  distribution 

MG  is  a  kind  of  small  power  distribution  and  consumption 
system,  which  overcomes  the  deficiencies  of  DG  in  terms  of 
intelligence  and  flexibility  [52].  MG  provides  power  energy  for 
users  by  modern  power  techniques,  such  as  fast  power  electronic 
switches,  advanced  converting  techniques,  efficient  new  energy 
sources  and  various  energy  storage  devices.  Many  kinds  of 
distributed  generators  and  energy  storage  devices,  including 
photovoltaic  (PV)  arrays,  wind  turbines  (WTs),  microturbines 
(MTs),  fuel  cells  (FCs),  and  energy  storage  batteries,  are  widely 
adopted  in  MG. 

MG  is  a  controllable  power  system  that  provides  thermal  and 
electricity  for  users  by  combining  the  power  generation  equip¬ 
ment,  energy  storage  devices,  control  equipment  and  load  [53], 
The  basic  structure  of  MG  has  been  presented  by  the  America 
Consortium  for  Electric  Reliability  Technology  Solution  (CERTS) 
[54],  MG  can  operate  in  both  grid-connected  mode  by  collaborat¬ 
ing  with  main  power  grid  and  independent  isolated  mode.  As  an 
effective  supplement  of  main  power  grid,  MG  has  attracted 
considerable  attention  in  recent  years  due  to  its  characteristics  of 
low  cost,  low  pollution,  high  reliability  and  easy  to  control. 

The  optimal  load  distribution  of  MG  is  to  achieve  multiple 
goals,  including  minimizing  cost,  minimizing  pollution  emissions, 
and  maximizing  reliability,  etc.,  while  satisfying  various  con¬ 
strains.  Fig.  1  shows  the  basic  structure  of  MG  optimal  load 
distribution. 

There  are  many  differences  between  the  MG  and  traditional 
large-scale  thermal  power  generation.  First,  power  generation 
of  various  distributed  generators  in  MG  usually  follows  the 
maximum  power-point  tracking  mode  [55,56],  and  cannot  be 
controlled  through  manual  scheduling.  The  output  characteristics 
of  these  distributed  generators  of  renewable  energy  power  gen¬ 
eration  are  sensitive  to  the  natural  conditions,  such  as  the 
intensity  of  sunlight  and  the  velocity  of  wind.  Second,  the  power 
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Fig.  1.  Basic  structure  of  MG  optimal  load  distribution. 


line  loss  of  MG  cannot  be  ignored  since  the  wiring  resistance 
caused  by  the  low  voltage  is  large.  Third,  the  equal  incremental 
principle  [57]  cannot  be  directly  used  for  MG,  since  the  power 
generation  characteristics  of  distributed  generators  in  MG  are 
different  from  traditional  thermal  power  units.  Also,  MG  should 
be  optimized  globally  and  dynamically  according  to  the  system 
operation,  since  the  power  distribution  of  MG  may  change  a  lot 
when  converting  between  grid-connected  mode  and  isolated 
mode  [58], 

Therefore,  traditional  optimal  load  distribution  models  of  the 
large-scale  power  systems  [59-61]  cannot  be  directly  applied  to 
MG  due  to  its  specific  characteristics.  New  optimal  load  distribu¬ 
tion  models  for  MG  need  to  be  developed. 

3.  Optimal  load  distribution  models  of  MG 

3.1.  Single-objective  models 

The  optimal  load  distribution  of  MG  is  actually  a  multi¬ 
objective  optimization  problem.  However,  in  some  existing 
research  work,  the  multi-objective  optimal  load  distribution 
models  of  MG  were  always  simplified  to  single-objective  ones  in 
order  to  facilitate  problem  analysis  and  solving.  Chen  and  Zhu  [58] 
established  a  single-objective  model  of  MG  optimal  load  distribu¬ 
tion.  The  objective  of  their  model  is  to  minimize  the  total  cost 
of  power  generation  by  optimizing  the  outputs  of  distributed 
generators  while  meeting  the  system  operation  conditions. 
The  objective  function  of  their  model  is  defined  as  follows: 

min  F=  min  £  [  2  Fi(Pi(t))+EbuyPbuy(t)-Ese„Psell(t)]  (1) 

t  =  1  i  =  1 

where  F  is  system  total  cost  of  power  generation,  T  is  the  total 
number  of  periods  in  scheduling  cycle,  t  is  the  number  of  period,  N 
is  the  total  number  of  power  generation  units  and  energy  storage 
devices  which  can  be  scheduled  within  the  system,  P,  represents 
the  actual  power  output  of  unit  i,F(P,)  denotes  the  operating  cost 
of  the  ith  power  generation  unit  or  energy  storage  device,  Pbuy  is 
the  electricity  bought  from  the  main  power  grid,  Pseu  is  the 
electricity  sold  to  the  main  power  grid,  Ebuy  is  the  price  of  bought 
electricity,  and  Eseu  is  the  price  of  sold  electricity. 

There  are  four  constraints  in  this  optimization  model,  including 
the  operation  constraints  of  power  generation  units,  charge- 
discharge  and  capacity  constraints  of  energy  storage  devices,  and 
power  balance  constraints.  The  total  operating  cost  in  this  model 


includes  fuel  cost,  maintenance  cost,  start-up  cost,  and  trading 
cost  and  income  with  main  power  grid,  but  does  not  include  the 
factor  of  fuel  consumption  rate.  Moreover,  the  objectives  of 
pollution  emissions,  system  reliability,  power  quality,  and  opera¬ 
tion  efficiency  were  not  considered  in  the  model.  Also,  the  time 
and  quantity  constraints  of  start-stop  were  not  included  as  well. 
Therefore,  the  simulation  results  of  this  model  may  not  be  so 
accurate. 

Hernandez-Aramburo  et  al.  [62]  pointed  out  that  the  most 
significant  difference  between  MG  and  large-scale  power  genera¬ 
tion  systems  is  the  presence  of  a  local  heat  demand,  which  adds 
another  dimension  to  the  optimization  problem  of  MG.  They 
proposed  an  optimal  load  distribution  model  that  considered  the 
Combined  Heat  and  Power  (CHP)  System.  The  aim  of  their  model 
was  to  reduce  the  fuel  consumption  rate  while  meeting  the  local 
energy  demand  (both  electrical  and  thermal).  However,  the  emis¬ 
sions  and  the  maintenance  costs  were  not  included  in  the 
objective  function  of  their  model.  Therefore,  it  is  difficult  to 
balance  the  relationship  between  cost  and  emissions  when  using 
their  model. 

Asanol  and  Bando.  [63]  and  Asano  et  al.  [64]  also  studied  the 
optimal  load  distribution  of  MG  considering  CHP  system.  The 
objective  function  of  the  mathematical  model  proposed  in  [63] 
includes  operating  cost,  gas  cost,  electricity  bought  cost  and  start¬ 
up  cost  of  gas  engine.  The  objective  function  proposed  in  [64] 
takes  more  costs  into  consideration,  including  gas  costs,  transac¬ 
tion  costs  with  main  power  grid,  the  contract  demand  cost  from 
the  main  power  grid,  operation  and  maintenance  costs,  start-up 
costs  and  the  original  investment  costs.  The  constraints  of  this 
model  include  power  supply  and  demand  balance,  charge-dis¬ 
charge  constraints  of  energy  storage  devices,  space  cooling  con¬ 
straints,  heat  gas  and  hot  water  demand  constraints.  However, 
both  the  two  models  [63,64]  were  simplified  to  single-objective 
optimization  problems.  Moreover,  capacity  constraints  of  energy 
storage  devices  and  operation  constraints  of  power  generation 
units  were  not  included. 

Zoka  et  al.  [65]  studied  the  optimal  operational  issues  of 
distributed  generators  and  energy  storage  devices  from  the  user 
perspective,  and  proposed  an  economic  evaluation  method  for 
independent  network  of  distributed  generator  resources.  They  also 
pointed  out  that  the  total  cost  of  MG  should  include  not  only  the 
installation  and  operating  costs,  but  also  self-construction  cost  and 
power  interruption  cost  which  are  the  measures  of  MG  reliability. 
The  total  cost  in  the  objective  function  is  composed  of  operating 
costs,  investment  cost,  and  the  revenue/costs  of  the  transactions 
with  main  power  grid.  The  operating  costs  consist  of  fuel  costs  and 
inspection  and  maintenance  costs.  The  constraints  in  the  model 
include  power  supply  and  demand  balance  constraint,  heat  supply 
and  demand  balance  constraint,  charge-discharge  balance  con¬ 
straint  of  energy  storage  devices,  generation  capacity  constraint, 
the  constraints  of  other  distributed  generators,  and  non-negative 
constraints.  The  advantages  of  this  model  are  twofold.  First,  its 
objective  function  integrates  the  concept  of  power  interruption 
cost,  which  indicates  the  reliability  of  MG  system.  Second,  the 
constitution  of  both  its  objective  function  and  constraints  is  more 
comprehensive.  However,  the  start-up  costs  and  the  line  loss  are 
ignored  in  the  objective  function.  Also,  energy  storage  capacity 
constraint  is  not  included  in  the  constraints. 

In  summary,  all  of  the  abovementioned  MG  optimal  load 
distribution  models  were  simplified  to  single-objective  optimiza¬ 
tion  problems.  However,  the  optimal  load  distribution  of  MG  is  a 
multi-objective  optimization  problem  in  which  we  should  not  only 
consider  the  objective  of  costs,  but  also  take  the  objectives  of 
emissions,  reliability,  power  line  loss,  power  quality  and  power 
generation  efficiency  into  account.  These  objectives  are  of  great 
importance  to  the  precise  and  effective  optimal  load  distribution 
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of  MG.  Therefore,  subjectively  ignoring  some  of  these  objectives  to 
simplify  the  multi-objective  model  into  single-objective  ones  is 
unreasonable,  which  can  directly  affect  the  accuracy  and  effec¬ 
tiveness  of  the  MG  optimal  load  distribution  results. 

3.2.  Multi-objective  models 

Mohamed  and  Koivo  [66-69]  studied  the  optimal  modeling 
and  online  management  of  MG,  and  established  a  series  of  multi¬ 
objective  nonlinear  mathematical  models  of  MG  optimal  load 
distribution,  in  which  the  minimizing  of  costs,  emissions  and  the 
trade-off  between  them  were  considered  at  the  same  time.  The 
operating  costs  of  the  multi-objective  optimization  model  estab¬ 
lished  in  [66]  are  composed  of  fuel  costs  and  operation  and 
maintenance  costs,  while  the  start-up  cost,  revenue  and  costs 
with  main  power  grid  are  not  included. 

The  total  emissions  of  pollutants  in  [66]  are  expressed  as 
follows: 

£(P,)=  Z  10-2(a/+M  +  riPf)  +  fiexpaiPi)  (2) 

i  =  1 

where  £(P,)  is  the  emissions  of  unit  i,  a, )q,Ci  and  A,-  are 
nonnegative  coefficients  describing  the  emission  characteristics 
of  the  ith  unit,  and  P;  is  the  real  output  of  the  ith  unit. 

The  constraints  in  this  model  are  inadequate,  which  only 
include  the  system  power  balance  constraint  and  distributed 
generator  operation  constraint.  Constraints  of  energy  storage 
devices  and  system  start-up  time  and  quantity  are  not  included. 
The  incomplete  model  will  lead  to  inaccurate  results. 

In  [67],  the  authors  analyzed  the  cost  characteristics  of  wind 
power  generation  systems,  photovoltaic  systems,  diesel  generators, 
fuel  cells,  micro  turbines,  and  batteries,  and  proposed  the  correspond¬ 
ing  cost  calculation  methods.  Then  a  multi-objective  optimization 
model  of  MG  was  established.  The  emission  objective  function  is  the 
same  as  Eq.  (2),  while  system  start-up  cost  is  added  to  the  operating 
cost.  However,  some  important  factors,  such  as  the  revenue  and  costs 
with  main  power  grid,  are  still  not  considered.  For  the  constraints,  the 
continuous  running/stop  time  constraint  and  start-stop  quantity 
constraint  are  added.  But  the  charge-discharge  constraint  and  capa¬ 
city  constraints  are  not  taken  into  account. 

The  optimal  load  distribution  model  established  in  [68]  is  much 
more  complete  than  the  previous  two  models  in  [66,67].  The 
emission  objective  function  is  the  same  as  Eq.  (2),  and  operating 
costs  include  fuel  cost,  operation  and  maintenance  cost,  start-up 
cost,  and  the  revenue  and  costs  with  main  power  grid.  The 
operating  costs  function  is  presented  as  follows: 

F(P.)  =  2  (CiFj  +  OMj+SCj  +  CPPj  —  ISPj)  (3) 

i  =  1 

where  F(Pf)  is  the  total  operating  cost  of  unit  i,  P,  represents  the 
actual  power  output  of  unit  i,  N  is  the  total  number  of  units  in  the 
MG  system,  Q  is  the  fuel  costs  of  unit  i,  F,  is  the  fuel  consumption 
rate  of  unit  i,  OM,  is  the  operation  and  maintenance  cost  of  unit  i, 
SCi  is  the  start-up  cost  of  unit  i,  CPP ,  is  the  cost  of  electricity 
purchased  from  the  main  power  grid,  and  ISP ,  is  the  income  of 
electricity  sales  to  the  main  power  grid. 

Constraints  of  the  model  in  [68]  are  the  same  as  those  in  [67], 
while  the  trading  between  MG  and  main  power  grid  is  considered 
more  comprehensively.  Two  different  conditions  are  considered  to 
model  the  purchased  and  sold  electricity. 

CPP  =  Cpx  max(Pt-Pi,0)  (4) 

lSP  =  Csx  max(P,  — Pj.,0)  (5) 

where  Cp  and  Csare  tariffs  of  the  purchased  and  sold  electricity 


respectively,  P,  is  the  generated  electricity,  and  PL  is  the  load 
demand. 

In  [69],  the  authors  added  the  emission  objective  measured  by 
emission  factor  to  the  total  objective  function.  Operating  cost  and 
constraints  are  the  same  as  those  in  [68],  This  model  also  has  some 
inadequacies.  For  example,  the  reliability,  power  line  loss,  power 
quality,  and  power  efficiency  are  not  considered  in  the  objective 
function,  and  the  charge-discharge  constraint  of  energy  storage 
devices  is  not  included  in  the  constraints  of  the  model. 

Wang  et  al.  [70]  studied  the  optimal  management  of  distrib¬ 
uted  generation,  and  presented  a  multi-objective  niche  evolution¬ 
ary  immune  algorithm  (MO-NEIA).  The  objective  function 
proposed  includes  electricity  quality  indicated  by  voltage  devia¬ 
tion,  operation  efficiency  of  generators  indicated  by  electricity 
factor,  and  power  line  loss.  The  constraints  include  variable 
constraints  (active  power,  reactive  power  and  voltage  constraints), 
and  power  flow  equation  constraints.  The  main  advantage  of  this 
model  is  that  it  takes  the  power  quality,  the  operation  efficiency  of 
generators,  and  the  power  line  loss  into  account,  and  the  quanti¬ 
tative  representations  of  these  factors  are  given.  These  are  rarely 
considered  in  other  models.  However,  some  costs  are  still  not 
considered  in  their  objective  function,  and  system  power  balance 
constraint,  charge-discharge  and  capacity  constraint  are  not 
considered  in  constraint  conditions. 

Though  multi-objective  optimization  models  of  MG  optimal 
load  distribution  have  been  established  in  the  aforementioned 
studies,  the  main  objectives  are  still  operating  cost  and  emissions. 
Many  important  objectives,  such  as  reliability,  power  line  loss, 
power  quality,  power  generation  efficiency,  etc.,  are  still  ignored. 
Meanwhile,  the  constraints  are  also  inadequate. 

Based  on  the  review  and  analysis  of  existing  studies  on  MG 
optimal  load  distribution,  we  find  that  many  single-  and  multi¬ 
objective  optimization  models  have  been  proposed.  Some  studies 
simplified  the  multi-objective  optimization  problem  to  single¬ 
objective  one.  Also  some  multi-objective  models  established  were 
still  inadequate  in  terms  of  objective  functions  and  constraints.  All  of 
these  problems  will  affect  the  precision  and  effectiveness  of  the 
optimal  load  distribution  results  of  MG.  Therefore,  a  more  compre¬ 
hensive  model  for  MG  optimal  load  distribution  model  is  needed. 


4.  A  comprehensive  model  of  MG  optimal  load  distribution 

4.3.  Objective  function  of  MG  optimal  load  distribution 

Actually,  the  optimal  load  distribution  of  MG  is  a  multi¬ 
objective  optimization  problem.  Though  the  emission  of  MG  is 
not  so  much  as  the  traditional  thermal  power  generation  since 
renewable  power  generations  are  widely  used  in  MG,  we  still 
cannot  ignore  the  factor  of  emission  in  the  objective  function.  In 
addition,  a  key  characteristic  of  MG  is  its  reliability  which  can  be 
measured  by  the  power  interruption  costs.  Also,  the  electricity 
generated  by  MG  is  low  voltage,  so  the  power  line  loss  cannot  be 
ignored.  Moreover,  since  transactions  exist  between  MG  and  the 
main  power  grid  when  power  generated  by  MG  is  not  equal  to  the 
load  demand,  the  revenue  and  cost  of  transactions  between  MG 
and  the  main  power  grid  should  be  included  in  the  total  operating 
cost.  Finally,  power  quality  and  operation  efficiency  of  MG  should 
also  be  considered.  Therefore,  the  objective  function  of  MG 
optimal  load  distribution  is  defined  as  follows: 

min  F(P,)  =  min  lFr(Pi),F2(Pi),F3(PilF4(Pi),F5(Pi),F6(Pi)}  (6) 

Fi(Pi)=  2  (CjFj  +  M,-  +  Sj  +  CEj  +  CP, - CSj) 


(7) 
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F2(Pi)  = 

UJ 

ii 

(8) 

N 

F3(P,)  = 

Z  ^ 

i=  1 

(9) 

N 

F4(Pf)  = 

Z  U 

i  =  1 

(10) 

N 

Fs(Pi)  = 

Z  Q, 

i  =  1 

01) 

N 

Fs(Pi)  = 

Z  Eff 

i=  1 

(12) 

where  F(Pj)  is  the  overall  objective  of  MG  system,  Fi(Pj)  is  the  total 
operating  costs,  F2(Pj)  is  the  total  emissions,  F3(P,)  is  the  interruption 
cost  which  indicates  the  reliability  of  MG,  F4(Pi)  is  the  total  power 
line  loss,  F5(P,)  is  the  cost  of  maintaining  an  acceptable  power 
quality,  andFg(Pi)  is  the  total  power  generation  efficiency.  In  each  of 
the  objective  function,  P,-  represents  the  actual  power  output  from 
unit  i  in  the  MG  system,  N  is  the  total  number  of  units  in  MG,  C;  is 
the  fuel  cost  of  unit  i,  F,-  is  the  fuel  consumption  rate  of  unit  i,  Mt  is 
the  maintenance  cost  of  unit  i,  S,-  is  the  start-up  cost  of  unit  i,  CEf  is 
the  emission  cost  of  unit  i,  CP,-  is  the  cost  of  electricity  purchased 
from  the  main  power  grid  when  the  load  demand  exceeds  the 
generated  electricity  from  MG,  CS;  is  the  income  of  electricity  sold  to 
the  main  power  grid  when  the  generated  electricity  from  MG 
exceeds  the  load  demand,  £,-  is  the  emission  of  unit  i,  /s’,  is  the  power 
interruption  cost  of  unit  i  which  indicates  its  reliability,  I,  is  the  line 
loss  of  unit  i,  Q,  is  the  cost  of  maintaining  an  acceptable  power 
quality  of  unit  i,  and  Eff,  is  the  operation  efficiency  of  unit  i.  For  the 
distributed  generators  in  MG,  their  operation  efficiency  “Eff  ’’  is  the 
efficiency  of  converting  the  input  energy  into  electrical  energy.  For 
the  energy  storage  devices  in  MG,  their  “Eff  means  the  power  output 
efficiency  of  the  stored  electrical  energy.  In  general,  the  operation 
efficiency  of  MG  refers  to  the  total  energy  input-output  efficiency. 

See  related  references  for  the  calculation  of  each  optimization 
objective. 

4.2.  Constraints  of  MG  optimal  load  distribution 

4.2.1.  Power  supply  and  demand  balance  constraints 

Power  system  balance  constraints.  The  MG  system  should  keep 
balance  between  power  supply  and  demand.  This  equality  con¬ 
straint  is  denoted  as  follows: 

N  Ns 

2  Pi+  ESj  +  PpUrchase  —  Psell  =  Pl+L  (13) 

i  =  1  i  - 1 

where  P;  is  the  actual  output  of  unit  i,  ESt  is  the  stored  energy  of 
energy  storage  unit  i,  N  andNs  are  the  total  number  of  generators 
and  energy  storage  devices  respectively  in  the  MG  system, Ppurdlase, 
is  electricity  purchased  from  the  main  power  grid,  Pse(/  is  the 
electricity  sold  to  the  main  power  grid,  PL  is  the  total  load,  and  Lis 
the  total  power  line  loss. 

Generation  capacity  constraints.  The  actual  output  of  each  unit 
should  satisfy  its  lower  and  upper  output  limits  as  follows: 

Pf"  <  Pi  <  P™51  (14) 

where  P™n  are  P™™  are  the  lower  and  upper  power  out  limits  of 
unit  i,  respectively. 

4.2.2.  Operation  constraints  of  energy  storage  devices 

Charge-discharge  constraints.  Energy  storage  devices  play  impor¬ 
tant  roles  in  MG,  their  operation  must  satisfy  the  charge-discharge 


constraints. 


<  EC,-  <  EC I™ 

(15) 

<  EDj  <  ED f  “ 

(16) 

where  EC;  and  ED,  are  charged  and  discharged  amount  of  energy 
storage  device  i,  respectively;  EC™", EC"™  and  ED^.ED™*  are  the 
lower  and  upper  limits  of  charge  and  discharge,  respectively. 

Capacity  constraint.  The  capacity  of  energy  storage  devices  must 
be  in  a  certain  interval  range. 

£S|™n  <  ES;  <  £SjIlax  (17) 

where  ES,  is  the  capacity  of  energy  storage  unit  i,  ESj11"1  and  ES™* 
are  the  corresponding  lower  and  upper  limits. 

4.2.3.  Power  generation  constraints  of  renewable  energy 

PV  power  constraint.  The  output  of  PV  power  is  significantly 
affected  by  natural  conditions,  such  as  the  temperature  and  the 
solar  radiation  at  a  certain  time  and  in  specific  region.  The  PV 
power  generated  cannot  exceed  the  power  capacity  under  stron¬ 
gest  nature  condition. 

Ppvj<PpvX  08) 

where  PPVi  is  the  power  generated  by  PV  unit  i,  and  P™x  is  its 
upper  limit. 

Wind  power  constraint.  Wind  velocity  has  a  great  impact  on 
wind  power  generation.  The  wind  power  generation  is  limited  as 
follows: 

PWT,i<Pw?  (19) 

where  PWT,i  is  the  power  generated  by  wind  power  unit  i,  and  P^ff 
is  its  upper  limit. 

4.2.4.  CHP  system  constraints 

CHP  power  constraint.  The  power  generated  from  CHP  system 
cannot  exceed  the  maximum  value,  which  is  expressed  as  follows: 

Pchpj  <  Penh  (20) 

where  PCHP,i  is  the  power  generated  from  CHP  unit  i,  and  P™pi  is 
its  upper  limit. 

Continuous  running/stop  time  constraint.  Once  the  units  are 
switched  on,  they  have  to  operate  continuously  for  a  certain 
minimum  time  before  being  switched  off  again.  On  the  other 
hand,  they  have  to  maintain  the  termination  status  for  a  certain 
time  before  being  started  again.  The  violation  of  these  constraints 
may  shorten  the  life  of  the  equipment  [68].  The  continuous 
running/stop  time  constraints  are  described  as  follows: 


zfon 
\lt- 1,; 

i-MUTiX0t_y-0t,i)Z  0 

(21) 

<f°/Lu 

—  MDT i)(()r-  i  j  —  0t  i)  >  0 

(22) 

where  T“"  ,  t/Tf[ , ,  is  the  on/off  time  of  unit  i,  MUT ,  and  AIDT,  are 
the  minimum  running  and  top  time  limits  of  unit  i  respectively, 
and  (9t_! ;  represents  the  on/off  status  of  unit  i  represented  by  1/0. 

Quantity  constraint  of  start-stop  units.  The  quantity  of  start  and 
stop  units  cannot  exceed  a  certain  number. 

N start  -stop  ^  hi  max  (23) 

where  Nstart_stop  is  the  number  of  start-stop  units,  and  Nmax  is  the 
upper  limit. 

4.2.5.  Other  constraints 

Reliability  constraint.  The  failure  of  MG  units  will  lead  to  the 
system  outage.  Taking  the  outage  probability  of  MG  into  account, 
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the  reliability  constraint  is  defined  as  follows: 

Prni  <  aPL  (24) 

where  Pm/  is  the  load  that  may  loss,  PL  is  the  total  load,  and  a  is  the 
probability  coefficient. 

Emissions  constraint.  The  distributed  generators  or  energy 
storage  devices  that  use  fossil  fuel  will  generate  a  certain  amount 
of  emissions.  Emissions  of  MG  must  be  limited  in  a  certain  interval 
range  according  to  related  environmental  protection  policies. 
Emissions  constraint  is  expressed  as  follows: 

I  Yu  <  yjnax  (25) 

!  =  1 

where  Yy  is  the  amount  of  pollutant/  generated  by  unit  i,  and  Yj1™ 
is  the  upper  limit  of  pollutant  j. 

Voltage  deviation  constraint.  The  voltage  deviation,  which  is 
typically  a  measure  of  power  quality,  of  distributed  generators  and 
energy  storage  devices  in  MG  cannot  exceed  an  upper  bound, 
which  is  described  as  follows: 

AY,-  <  AVmax  (26) 

where  AY;  is  the  voltage  deviation  of  unit  i,  and  AYmax  is  the 
upper  limit. 

User  defined  constraints.  In  addition  to  the  aforementioned 
various  constraints.  We  should  also  note  that  emissions  and 
efficiency  in  different  equations  may  result  in  dependent  rows  in 
the  array  to  solve.  Therefore,  in  practical  applications,  this  kind  of 
user  defined  constraints  should  also  be  considered. 


5.  Conclusion 

The  optimal  load  distribution  of  MG  is  an  important  aspect  to 
achieve  its  optimal  operation  and  efficient  application  in  the  smart 
grid  environment,  and  the  establishment  of  optimization  mathe¬ 
matical  models  is  the  first  step  of  MG  optimal  load  distribution. 
We  first  review  some  existing  optimal  load  distribution  models  of 
MG,  and  pointed  out  the  shortcomings  of  these  models.  Then  we 
propose  a  general  and  comprehensive  model.  It  is  of  great 
importance  to  improve  the  accuracy  and  effectiveness  of  MG 
optimal  load  distribution  models,  thereby  promoting  the  study 
of  MG  economics  and  achieving  the  intelligent  power  generation 
and  distribution  of  smart  grid. 

We  should  note  that  the  selection  of  an  appropriate  MG 
optimal  load  distribution  model  in  practical  smart  grid  applica¬ 
tions  is  important.  The  model  established  in  this  paper  is  a  general 
model,  and  we  should  consider  the  specific  characteristics  while 
using  the  model.  In  addition,  it  is  important  to  explore  more 
effective  and  efficient  algorithms  to  solve  the  optimal  load 
distribution  problems  of  MG  after  the  model  was  established. 
Some  of  the  traditional  optimization  methods  [71-73]  may  not  be 
suitable  for  MG  optimal  load  distribution,  which  is  a  nonlinear, 
multi-objective  and  multi-constraint  optimization  problem.  Some 
intelligent  algorithms,  such  as  genetic  algorithms,  bee  colony 
optimization,  and  particle  swarm  optimization  algorithm,  are 
promising  methods  to  solve  the  MG  optimal  load  distribution 
problems  [74-78], 
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